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We search for planar deviations of statistical isotropy in the Wilkinson Microwave Anisotropy 
Probe (WMAP) data by applying a recently introduced angular-planar statistics both to full-sky 
and to masked temperature maps, including in our analysis the effect of the residual foreground 
contamination and systematics in the foreground removing process as sources of error. We confirm 
earlier findings that full-sky maps exhibit anomalies at the planar (I) and angular (£) scales (I, £) = 
(2, 5), (4, 7), and (6,8), which seem to be due to unremoved foregrounds since this features are 
present in the full-sky map but not in the masked maps. On the other hand, our test detects 
slightly anomalous results at the scales {1,1) = (10,8) and (2,9) in the masked maps but not in the 
full-sky one, indicating that the foreground cleaning procedure (used to generate the full-sky map) 
could not only be creating false anomalies but also hiding existing ones. We also find a significant 
trace of an anomaly in the full-sky map at the scale (1,1) = (10,5), which is still present when we 
consider galactic cuts of 18.3% and 28.4%. As regards the quadrupole (£ = 2), we find a coherent 
over-modulation over the whole celestial sphere, for all full-sky and cut-sky maps. Overall, our 
results seem to indicate that current CMB maps derived from WMAP data do not show significant 
signs of anisotropies, as measured by our angular-planar estimator. However, we have detected a 
curious coherence of planar modulations at angular scales of the order of the galaxy's plane, which 
may be an indication of residual contaminations in the full- and cut-sky maps. 

PACS numbers: 98.80.-k, 98.70.Vc, 98.80.Es 



I. INTRODUCTION 

The statistical concordance model of the universe has now been established with unprecedent accuracy by the 
Wilkinson Microwave Anisotropy Probe (WMAP). The five year dataset release of the WMAP team PQH] has shown 
that - to a very good extent - the temperature field of the cosmic microwave background radiation (CMB) obeys a 
Gaussian statistic with zero mean and a variance which is the same in every direction in the sky. At the same time, 
several different teams (3j SI El EJ [Tl [HI El EH HU H2 US] have reported the detection of statistical peculiarities in this 
field, mainly in the largest cosmological scales, which range from slightly to extremely unlikely within the framework 
of a Gaussian and statistically isotropic universe. Among these anomalies, the most conspicuous are the lack of power 
in the low-£ multipole sector and in the 2-point correlation function P3IEL the alignment of the quadrupole £ — 2 and 
the octopole I = 3 E], and the so-called north-south asymmetry |U [HI El EB H2] [131 ESI [TBI E3 HH1 US]- Several 
attempts to explain away these statistical features in terms of known sources of microwave radiation or peculiarities 
of the instrument have been proposed, but have not yet produced a compelling explanation. 

If we put aside the hypotheses that these anomalies may possibly be due to some residual galactic contamination |20l 
I2T1 122] [23] or even to a systematic data analysis effect (24j 125] [26] . we inevitably end up reverse-engineering physical 
and/or astrophysical phenomena in order to explain these effects. Although the former explanations are important 
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and deserve attention on their own, the possibility that a new mechanism may explain some of the unknown CMB 
features is still far from being ruled out. In such a case, we have essentially two different approaches: a theoretical 
(bottom-up) and a phenomenological (top-down) one. 

In the bottom-up approach, models for the evolution of the universe are formulated through physical principles in 
such a way as to account for deviations of Gaussianity and/or statistical isotropy (SI). These models usually invoke the 
existence of non-trivial cosmic topologies (23 EH ES] , primordial magnetic fields (12 ED |32] |33l [34] , local (3l | l36 j l37l 138] 
and global [39] gQ] |HJ 51] manifestations of anisotropy and even exotic models of inflation g3] 23J g5] 26] g7] 135] . 
only to mention a few of the possibilities. 

In the phenomenological top-down approach, we design statistical tools to test the robustness of temperature maps 
against the hypotheses of Gaussianity and SI. These approaches include statistics of multipole alignments using 
Maxwell's multipole vectors [6] [39] |50] [51] , constructs of a posteriori statistics [8j and functional modifications of the 
two-point correlation function }52j [53] 1531 [55] . 

One interesting feature of the top-down approach is that it can shed some light on the source of the reported 
anomalies without the need to invoke any specific cosmological model, but rather acting as a further guide for theorists. 
However, due to the intrinsic generality of this approach, it is sometimes unclear how to select the important degrees 
of freedom in order to construct meaningful statistical estimators. This is in fact a natural consequence of model- 
independent approaches, which at some stage force us to rely on our theoretical prejudices about the statistical nature 
of the CMB in order to construct estimators of non-Gaussianity and/or statistical anisotropy. 

We can instead follow another route and let the construction of our statistics be guided primarily by what the 
observations tell us, rather than by what we think they should look like. Here we are primarily interested in answering: 
how to detect planar signatures in CMB maps? One possibility was recently proposed in [56], where the two-point 
correlation function is decomposed in such a way as to account for the presence of symmetry planes in temperature 
maps. Given that CMB experiments are confined to collecting data from inside our galactic plane, it is conceivable 
that we may find planar signatures in the data if foregrounds were improperly removed. Besides being observationally 
motivated, this approach leads naturally to both an unbiased estimator and a rotationally invariant test of statistical 
anisotropy, therefore alleviating the difficulties mentioned above. 

In this work we employ the angular-planar power spectrum analysis introduced in |56] to search for planar 
signatures both in partial and full sky-coverage temperature maps, allowing for possible effects due to residual 
foreground contaminations. Our analysis may therefore help to shed some light in the origin of the reported anomalies 
and test the robustness of cleaned maps due to foreground contamination in a model-independent way. 

This work is organized as follows: we start ^jn]by reviewing the construction of the two-point correlation function, 
and we show how it can include planar deviations of SI. We then show that a simple chi-square analysis can be 
designed to detect such planarities in any theoretical cosmological model. We end that Section with a simple example 
of planarity which could be easily detected through our tests. In pil| wc discuss how the probabilities of measuring 
particular values of cosmological observables are usually calculated, and how this calculations can be extended to 
include the uncertainty inherent in this measurement. In Section §IV| we explain how to include the effect of foreground 
contamination in our analyses, and present our results for full-sky temperature maps. After presenting the results of 



the analyses for masked maps in { IV B we conclude and give some perspectives of further developments. 



II. ISOTROPIC AND ANISOTROPIC TWO-POINT CORRELATION FUNCTIONS 

The starting point of our analysis is the temperature fluctuation field AT(n). According to standard inflationary 
models, this is a random field of which we have access to one single realization. If we decompose this field in spherical 
harmonics: 

AT(n) = ^2 a-imYim{n) , (1) 

its randomness can be attributed to the multipolar coefficients ae m , and the statistics of the field can be characterized 
by constructing the statistical moments of these coefficients. For simplicity, we suppose that this field is Gaussian 
with zero mean (AT(n)) = 0, in which case the two-point correlation function, 

C(h u n 2 ) = (AT(n 1 )AT(n 2 )) = £ £ K^a^J^^fni)^^^) , (2) 

is the only relevant statistical moment. Equivalently, the covariance matrix (af t m.i ct l 2TO2 ) encloses all the statistical 
information of the multipolar coefficients in Eq. ([l]). Notice that the hypothesis of Gaussianity makes no restriction 
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whatsoever upon the rotational properties of the two-point function C(hi,n2)- It is therefore perfectly legitimate to 
restrict our considerations to the Gaussian case when looking for deviations of SI, noting that these two questions 
should be clearly stated and, as much as possible, decoupled. 



A. Isotropic correlation function 

The statistical counterpart of the concordance ACDM model predicts a universe which, in its linear regime, is not 
only Gaussian but also statistically isotropic. This means that the temperature fluctuation field (and its variance) 
has no privileged directions, and that the variance depends only on the angular separation of two given points in the 
CMB sky: 

C(ni • na) = ^ 2 ^^C l P t {fi 1 ■ h 2 ) . (3) 

e 

In harmonic space, rotational invariance means that the correlation function is completely diagonal: 

(ae imi a*i 2m2 ) = Ct, 1 Se 1 e 2 5 mim2 , (4) 

with its eigenvalues determined by the theoretical angular power spectrum C? 1 . Together, Gaussianity and SI comprise 
a very severe set of restrictions to the temperature field Eq. ([TJ , and allows no expressions other than ([3} and |4]) to 
specify its properties. 



B. Planar correlation function 

Deviations of SI presuppose functional deviations of the two-point function from Eq. Despite the infinite 

possibilities that are open when one relaxes the requirement that the two-point function depends only on the 
separation angle cos 8 = hi ■ fi2, we should still remember that the unit vectors hi and h 2 share a common origin at 
the center of the CMB sphere (in other words, the center of our surface of last scattering.) One of the possibilities 
that do not assume isotropy is to consider C = C '(hi, hi) |52j, but in this case we would not actually be measuring 
correlations of two points, but rather the self-correlation of all individual points in the sky (i.e., the coincidence limit 
of the two-point correlation function.) Another possibility is to consider C = C(hi,h2) as given by expression 
with the two vectors completely independent of each other |53j. Unfortunately this function is too general, and leads 
to rampant arbitrariness in the construction of statistical estimators of anisotropies, as we mentioned in £JT| 

A third possibility was pointed in [56 : for any pair of unit vectors in the CMB sphere, not only their angular 
separation is well-defined, but also the plane where they stay. That is, since these vectors are restricted to have a 
common origin, both their scalar and cross products are uniquely defined. This suggests that we incorporate this 
vectorial dependence into the two-point correlation function: 

C(hi,h 2 ) = C(hi x h 2 ) , (5) 

which can now account for both angular and planar modulations of the temperature field. Note that now vectors with 
the same angular separation are not necessarily correlated in the same way. This feature is what makes the functional 
form above suitable for detecting planar signatures in CMB sky maps, no matter if they have an astrophysical origin 
(as in, e.g., the galactic plane) or if they are artificially imprinted on the maps (e.g., poor masks, inhomogeneous 
observation function, etc.) 

In practice, just as in the case of the usual (isotropic) C/s, it is easier to work with ([5} in harmonic space. If we 
introduce the definition n = hi x h 2 = {sin d, 8, <fi}, where $ = arccos hi ■ h 2 , we obtain the following decomposition: 

C(n) = VV^^Cf l F £ (cos^y im (n), I e 2N , (6) 

where the restriction of the sum above to even values of I is a consequence of the reciprocity relation C(hi,h 2 ) — 
C(h 2 ,hi). The decomposition above is equivalent to considering modulations of the function Cg(n) in the sphere, 
where the usual angular power spectrum Ci is now given by the monopole C°°. In fact, we could have started with 



4 



this modulations rather than (|6f , but in such a case we would perhaps neglect the geometrical interpretation of the 
normal vector n, which is crucial in our analysis. 

By equating Eqs. (JgJ and it is possible to show [56] that the coefficients C l e m are directly related to the 
temperature multipolar coefficients ag, m . We have that: 



l\ ,mi l 2 ,m 2 



f = 2n E E k^^) ( l m m \ m 2 2 J . ( 7 ) 



where 



4^ ee ]T(-irA, im A, 2m r<-cosi?)P,(co S i?)e^ ( J * 

m \ 



4 ^2 
TO —TO 



(8) 



with A^ m a set of real coefficients resulting from the d integration and which are zero unless I + m = even (see 
the Appendix of [56] for more details). These expressions show, first, that the angular-planar multipolar coefficients 
C^ m can be calculated from first principles in any theoretical model which predicts a specific form for the covariance 
matrix {ae 1 mi a e 2 m 2 ) ^ or instance, models of anisotropic inflation |39 | l4"0 | I4"T ] |4"2]. And second, given any observed 
temperature map, we can estimate the values of all C l e m, s by computing the sum in Eq. (|7j with the actual coefficients 
a>lm - see the next Subsection. Note also that the coefficients Q are not restricted to temperature maps, and can be 
equally applied in the analysis of CMB polarization, large-scale structure maps or in fact any map on S . 



C. Statistical estimators of anisotropy 

The multipolar angular-planar coefficients ([7| were defined in terms of an ensemble average of the temperature 
multipolar coefficients ai m . Of course, we only observe one universe, which means that the cosmic variance in the 
determination of the coefficients C™ is a severe restriction that we have to live with. This means that, if we want to 
evaluate the statistical properties of these coefficients, like mean values or variances, we have to build functions which 
can at best estimate these quantities. One obvious possibility is the following unbiased estimator: 

CjT "> 2 7 ^^/2^TI £ £ a lxMm ^ ( ^ * £ ) i& , (9) 

ii,mi l 2 ,m 2 

which can be uniquely calculated for a given temperature map. 

Following this prescription, we can now estimate the statistical moments of (|9}. In particular, its mean value was 
already calculated, and is given by \/2l + 1 times the right-hand side of ^7§. Notice that, under the hypothesis of SI, 
it follows from Eqs. and Q that [56] : 

(C l " l )si = ClSioSmO . (10) 

Conversely, if the only non-zero C l e m, s are given by I = m = 0, then C®° — Ct, and we conclude that SI is achieved if 
and only if the coefficients C e m, s are of the form (10 1. From now on, we will restrict our analysis to / ^ 0, since this 
is the first non-trivial case of planarity. 



D. A x test of statistical isotropy 

We have determined above, in a statistically isotropic universe, that all C l e m 's with I ^ are random variables with 
zero mean. We can now ask, with a set of C^ m 's at hand, how do we check whether they respect statistical isotropy - 
and in case they do not, which type of anisotropy they correspond to? 



If we have a set of zero-mean random variables as in (10 1, an obvious choice would be to construct the associated 



chi-square test in the same way as the usual Cgs are constructed as a chi-square fit from the a£ m 's. However, there 



is nothing which prevents us from applying this test to models where ( 10 1 does not hold - all that is needed is to 
subtract from the C^ m 's its expectation value given by a fiducial model. Hence we can define the following (reduced) 
chi-square test: 

-i ^ \nlm lnlrch\ I 

r v 2y - 1 v 1 1 K 1 ;| mi 
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where (C l e m ) and (a 1 ™) 2 are respectively the expectation value and the variance of the C l e m, s in some particular model. 
Notice that, since we are summing over m in the definition of the (xt)e, our test is rotationally invariant. This is a 
crucial property, since we are pursuing a "blind test" of anisotropy where we do not know what type of directionality 
we are looking for. 



Now that we have defined a general test of anisotropy, we need a model to be tested. As it turns out, the most 
successful model we have is the concordance ACDM model, for which SI holds. Under the hypothesis of SI, not only 
(C 1 " 1 ) can be easily calculated, but also the variance (a 1 ™) 2 have a simple, m- independent, expression: 



(°i m )si = ( a e) 2 = ^ ^ Ct- c t- 



• 2 E 



(12) 



For this model, it is clear from Eqs. (10), (11 1, and (12| that we will have {{x 2 )\) = 1 as long as I ^ 0, so we define 
the angular-planar x 2 measure of the deviations of any map from SI to be: 



{X% = (xlYe 1 • 



(13) 



If the test x 2 ls significantly positive or negative (allowing for some variance which can be different for different maps) 
this will be an indication of deviations from SI. 



E. An example of planarity 



Before we move on to apply ( 13 1 to actual CMB data, it is important to have some intuition on the type of planarity 
which could be revealed by this test. We can ask, for example, what pattern of modulations would a Gaussian and 
SI temperature map AT(n) inherit when we apply to it a window function W(h). Depending on the shape of this 
function, the resulting map 



AT (n) = AT(n)W(n) 



(14) 

will contain planar modulations which our test ( |l3"| can easily detect. Notice, however, that this procedure is not the 
same as 'cutting' or making a 'hole' in the sphere S 2 in order to neglect the contaminations from our galaxy with 
the help of a mask (as we will do in SIVBl. The effect of the window function above is conceptually simpler, and 



corresponds to a region with enhanced/suppressed fluctuations, which naturally induces strong anisotropic correlations 
on AT. In particular, the resulting angular power spectrum Ct will not be trivially related to the isotropic CVs, but 
rather it will be given by a convolution with the latter which couples different angular scales through some kernel 
[361 157] . Consequently, even the angular spectrum Cg would hint strongly towards anisotropy, although in this case 
no further information of planarity could be extracted from it. 




Figure 1: Temperature map with a planar signature around the equator. See the text for more details. 



Having these distinctions in mind, let us take as an example a simple window function which gives zero weight for 
fluctuations around an equatorial strip of fixed width A9, i.e.: 

0e§[7r±A0] 

otherwise . 
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This window function introduces a disk around the equator where all points have the same temperature, as can be 
seen in Fig. (JTJ . If we further decompose the function ( |l5| in spherical harmonics and use (jlj , we arrive at (36] : 



E 

I'm 1 



TfZ m „ 



(16) 



with K\™ given by 



K, 



I'm 



^2w LM / d 2 nY LM (n)Y e > ml (h)Yl 

L,M J 



■in) 



and where wlm are the coefficients resulting from the decomposition of (15 I into spherical harmonics 



Although, in principle, expression (161 could be used to explicitly evaluate both (12 1 and (13 1, the resulting 



expression is not clarifying. In practice it is much easier to obtain the same result by simulating these maps and 
applying the x 2 test to each one of them. Nonetheless, the relevance of Eq. (161 lies in the linear relation that 
is established between ag m and a£ m , which shows that the resulting temperature field AT will still be Gaussian, 
provided that AT is Gaussian. In other words, in this example we have a Gaussian but anisotropic field (AT), and 
therefore we will be performing a test of SI only, which is completely decoupled from the issue of Gaussianity. 



Following the above prescription, we have applied the estimator (13) to 2 x 10 maps of the form (16 1, with the 



window function chosen to cover an equatorial strip 50° wide; the result is shown in Fig. fl2|. The first thing we notice 
in this figure is that small planar scales (i.e., a oc l^ 1 <c 1, where a is the angle formed by the normal of the planes 
that lie inside the disk and the z axis) have a roughly constant planar modulation for all the angular scales (reflected 
in i) we tested. This is indeed what we would expect, for in this case we are probing the temperatures at points lying 
in circles which are nearly parallel to the equator where, by construction, the temperature is constant. 

As we examine larger and larger planar scales (smaller I's), we see that the modulations become constant only 
below a certain angular threshold, corresponding to those vectors that still live inside the constant temperature strip. 
Hence, for the planar scales I = 2, 4 and 8, the angular correlations level out at angular scales of approximately I ~ 6, 
10 and 12, respectively. Notice that the plateau of the % 2 for each at i ~ 10 - 12, becomes smaller as I grows. This 
reflects the fact that the 50° disk induces correlations on large planar scales. These correlations result in a stronger 
signal in x 2 f° r the small angular scales (£ ^ 1) simply because for smaller angular separations there are more pairs 
of points that contribute coherently to % 2 within the disk. 



0.3 



0.2- 



0.1 



0.0- 




10 



12 



14 



Figure 2: The angular-planar \ 2 test averaged over 2 x 10 3 maps with a planar modulation of ±25° around the equator. The 
figure shows ((x 2 )«) as a function of £ for some particular scales in the range I £ [2, 20]. 



For a statistically isotropic map we would expect the (x 2 )^ to be randomly distributed around zero. But in the 
example above we see that this does not happen: there is a coherence of the (x 2 ) l £ over many I's and ^'s, denoting 
the preferred directions/planes. Hence, the main lesson from this example is that angular-planar features pointing 
towards statistical anisotropy can sometimes be individual anomalous values of some (x 2 )^'s, which would correspond 
to angular-planar anomalies at particular angular scales i and planar scales I. However, we can also have scenarios 
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where the individual (x 2 )^'s are n °t particularly anomalous, but the coherence of the test over a range of £'s and I's 
is an indication of some preferred regions or directions. The first situation [an individual anomalous (x 2 )^] points 
towards a sharp anisotropic feature, whereas the second scenario [our present example, of a coherent feature in all the 
(x 2 )^'s] points towards a wider region of anomalous correlations. 

To conclude this example, we mention that although it was artificially built to illustrate the use of the angular- 
planar estimator, we will see in §IV B| that similar modulations are found in cleaned CMB maps, even after applying 
the more realistic masks KQ75 and KQ85. 



III. CALCULATING PROBABILITIES 



Given that our main goal is to determine how "typical" our universe is according to the test (13 1, we will now 
revisit the question of the probabilities for such tests in cosmology. 

Let us suppose that we want to calculate the probability of detecting an observable inside a specific range of 
values, according to a given cosmological model. For instance, we may genuinely ask what is the probability that in 
a "randomly selected" Gaussian and SI universe we would measure one of the (x 2 )^ with a value equal or smaller 
(greater) than the value which is actually measured. Often, this question is analyzed as follows: suppose that the 
probability density function (pdf) of a given observable x in this particular model/theory is known and given by Vth- 
Let us suppose further that the observation of this quantity gave us a value xq. Then, the number: 

P< = I ° V th {x) dx (17) 

J —oo 

denotes the probability of detecting a value x equal to or smaller than xq, according to this model. The probability of 
having x > xq is then given by P > = 1 — P<, However, these probabilities assume that we have measured the value 
xq with infinite precision, so they say nothing about an important question that we must deal with: typically, the 
measurement of xo has itself an uncertainty, which should be folded into the final probability that the observations 
match the theoretical expectations. 

Since no measurement (of the CMB or any other physical observable) will give us a result with infinite precision, 
our measurements should also be regarded as random events. Therefore, in a more rigorous approach, we would 
have to consider xq itself as a random variable. In the case of CMB, however, this would be only part of the whole 
picture, since the randomness of the measurements of xo should also be related to the way that this data is reduced 
to its final form. This happens because different map cleaning procedures will lead to different values for xq. This 
difference induces a variance in the data which reflects the remaining foreground contamination of the temperature 
maps. We will come back to this point later, after we show how to calculate probabilities when both x and xq are 
random variables. 



A. Difference of random variables 



Suppose that x and y are two random variables, uniformly distributed over the range [0,1]. What is the pdf of 
the new variable z — x — yl Intuitively, we would expect this pdf to be zero for z = ±1, since the pairs of events 
(x, y) = (1, 0) and (0, 1) are the least frequent ones. By the same reasoning, we would also expect a peak around zero. 
The correct answer is a triangular pdf for z: 

V(z) = l-\z\ , 

which can be calculated by convolving the pdf 's of x and y. This result can be generalized to include any type of 
continuous random variables [58]. This is precisely what we need, and if we define: 

x = (X 2 ) l i(th) > V = (x 2 )^ ( obs) . z = x-y, 
then the pdf of z will be formally given by 

/oo poo 
V ob s(y)V th (z + y) dy = V ohn (x - z)V th (x) dx . (18) 
-oo J — oo 

Notice that this probability density will be automatically normalized if both Pobs an d Pth are - 
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According to our definitions, the probability that in a randomly selected ACDM universe we detect a value x < y 
is translated into the probability of having z < 0, and this is simply the area under V(z) for z £ (—00, 0]: 

P<= f V(z) dz . (19) 
J —00 

Moreover, in the ideal case where observations are made with infinite precision, Vobsiu) — Hu ~ x o) and we recover 



171) . The reader must be careful, though, not to think of (171 as a lower bound to (19 1, since the distributions V b S 



and Vth are not necessarily symmetric about their mean values. In other words, a large distance from y to the most 
probable theoretical value x would not, by itself, constitute sufficient grounds to claim that this measurement of y is 



'unusual', since its dispersion can be wide enough to render it 'usual' according to (19 1. 

Finally, there remains the question of how to obtain the probability densities T-'obs and 7-th- These functions can be 
computed numerically, provided that the number of realizations of the random variables y and x is large enough, since 
in this case the histograms for these variables can be considered as piecewise constant functions which approximate 
the real pdf 's. For the case of the ('theory') variable x defined above we have run 2 x 10 4 Monte Carlo simulations 
of Gaussian and statistically isotropic CMB maps using the ACDM best-fit GYs provided by the WMAP team [59 . 
With these maps we have then constructed 2 x 10 4 realizations of the variable x. This procedure was also carried 
using masked C/s (corresponding to masks KQ85 and KQ75) to construct 2 x 10 4 masked versions of the variable x. 

The simulation of the variable y, related to the observations, is less trivial, and is intrinsically related to the way 
we estimate contamination from residual foregrounds. We will now detail this procedure. 



IV. ESTIMATING THE RESIDUAL FOREGROUNDS 



As mentioned before, cosmic variance is an unavoidable limitation which can be estimated by, e.g, running many 
different realizations of our theoretical random variable x. However, cosmic variance is far from being the sole source 
of uncertainties in CMB experiments. As is well-known, not only instrumental noise, but systematic errors (e.g., in the 
map- making process), the inhomogeneous scanning of the sky (i.e., the exposure function of the probe), or unremoved 
foreground emissions (even after applying a cut-sky mask) could corrupt - at distinct levels - the CMB data. 

Foreground contamination, on the other hand, may have several different sources, and consequently several different 
ways of being included, many of which would certainly go beyond the scope of the present work. However, since 
different teams apply distinct procedures on the raw data in order to produce a final map, we will make the hypothesis 
that maps cleaned by different teams represent - to a good extent - "independent" CMB maps. We will use this 
hypothesis in order to estimate residual foreground contaminations by comparing these different foreground-cleaned 
maps. 

As a matter of fact, the WMAP science team has made substantial efforts to improve the data products by 
minimizing the contaminating effects caused by diffuse galactic foregrounds, astrophysical point-sources, artifacts from 
the instruments and measurement process, and systematic errors |60| IST] . As a result, multi-frequency foreground- 
cleaned full-sky CMB maps were produced, named Internal Linear Combination maps, corresponding to three and five 
year WMAP data [T][62]. Therefore, in order to account for the mentioned randomness, systematic, and contaminating 
effects of the CMB data, we include in our analyses several full-sky foreground-cleaned CMB maps, listed in Table [I] 
which were produced using the three and five year WMAP data. 



Full-sky maps 


References 


Hinshaw et. al. 


EH] 


de Oliveira-Costa et. 


%i. m] 


Kim et. al. 


m 


Park et. al. 


m 


Delabrouille et. al. 


m 



Table I: Full-sky foreground cleaned CMB maps from WMAP data used in our analysis to estimate a lm . Note that the 
reference |63j includes the analysis of maps from the three and five years WMAP releases. 

The prescription we adopt to determine the distribution of the observational variable y can be equally applied 
either to partial or full-sky CMB maps, and is achieved as follows: we simulate Gaussian random a^ m 's in such a way 
that their central values are given by the five year ILC5 data [TJ [62], and with a variance which is estimated from 
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the sample standard deviation of all the maps listed in Table |TJ So, for example, suppose we have n different full-sky 
temperature maps at hand and we want to estimate the randomness inherent in the determination of, let's say, 032- 
Therefore, we take: 



M{af 2 c \a Z2 ) 



032 



(20) 



with: 



032 



\ 



1 



(a\ 2 - a 32 ) 2 



and 



032 



i=l 



- t 32 > 



(21) 



where AT(n, a) represents a Gaussian distribution with mean \i and standard deviation a. 

Our procedure may be justified by the following facts: first, the ILC5 data is undoubtly the most accurate CMB 
measurement we currently have at our disposal, and for this reason it can be considered as a good approximation to 
the "real" temperature map. Secondly, as we already mentioned, by comparing different CMB maps we should have 
an estimate of the residual contamination which were not properly removed, and that can possibly be the source of 
anomalies. Note that if this contamination is indeed weak, then the sample variance above will be small, and our 



procedure will reduce to the standard way of calculating probabilities (see the discussion of [III I. Finally, the use of 



a Gaussian in (20 1 was dictated not only by simplicity, but also by the fact that the propagation of uncertainties in 



physical experiments is usually assumed to follow a normal distribution. 



A. Full-sky maps 

Following this procedure, we have used the full-sky maps shown in Table |l} to construct 10 4 Gaussian random ae m 's, 
which were then used to calculate 10 4 realizations of y = (x 2 )f( bs)- With those variables we constructed histograms 
which, together with the histograms for the (full-sky) variable x, were used to calculate the final probability ( [T9| . We 
have restricted our analysis to the range of values (£, I) S [2,10], since the low multipolar sector (i.e., large angular 
scales) is where most of the anomalies were reported. The resulting histograms and pdf's are shown in Fig. (|3j, and 
the final probabilities we obtained are show in Table (p5|. 
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10 
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81.1% 


73.7% 


54.1% 


6.1% 


80.7% 


46.4% 


36.9% 


47.5% 


81.8% 


4 


74.0% 


72.6% 


55.0% 


39.6% 


74.2% 


93.2% 


51.6% 


55.4% 


56.3% 


6 


78.1% 


80.7% 


69.3% 


52.3%, 


33.6% 


80.0% 


95.0% 


50.3% 


82.2% 


8 


63.5% 


87.7% 


18.8% 


51.5% 


21.4% 


66.4% 


31.6% 


27.5% 


82.3% 


10 


67.9% 


50.0% 


61.0% 


8.7% 


37.7% 


59.5% 


36.6% 


29.2% 


35.7% 



Table II: Final probabilities of obtaining, in a random ACDM universe, a chi-square value smaller or equal to (x 2 )^(obs)i as 
given by full-sky temperature maps. 

It is interesting to compare the values in Table |n} with the more naive and independent 'error-bars' estimate made 
in [56]. As anticipated there, the angular quadrupole I = 2 does not seem anomalous under the tests we considered. 
The only curiosity lies in the fact that, for the scales we probed, the quantity (x 2 ) 2 has a probability which is always 
greater than 50%, meaning that the quadrupole 1 — 2 has a consistently positive planar modulation. 

By far, the most unlikely individual values in Table (JllJ) are in the sectors (l,£) given by (2,5), (10,5), (4,7) and 
(6,8), and are all below a relative chance of 10% of either being too negative [(2,5), (10,5)] or too positive [(4,7), 
(6,8)]. Although these modulations do not extend to a significant range of I values, it is possible that these sectors 
are contaminated by foreground residuals. If this is the case, then these anomalies should not be present in masked 
temperature maps. Conversely, if these are the only anomalies present in CMB data, then the next cut-sky maps 
analyses should not reveal other extreme probability values (i.e., lower than 10% or higher than 90%). We will analyze 
these questions in more detail in the next Section. 

But before we turn to the masked maps, it is important to realize that not only the individual values of (x 2 )^ 
are relevant: their coherence over a range of angular or planar momenta also carries interesting information. So, 
for example, a set of (x 2 )^'s which are all individually within the cosmic variance bounds, but which are all positive 
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Figure 3: Full-sky maps pdf 's for (x 



^(th) 



(blue histograms) 



«(th) ~~ 

(x^iobs) (solid re d line). We show only a few representative figures since the remaining ones are qualitatively the same. The 



final probabilities are shown in Table and correspond to the area under the solid curve from 
normalized to 1. 



f(obs) (P ur pl e histograms) and for the difference (\ 
:e the remaining ones are quali 

oo to 0. All pdf's are 



(or negative) can be an indication of an excess (or lack) of modulation. In fact, as we have shown in SHE if one 



suppresses a Gaussian and isotropic map only over a thick disk between fixed latitudes, the effect is to boost/depress 
the modulations compared with the isotropic case, over a large range of Vs and £'s. The I's and £'s which are most 
affected are the ones corresponding to the angular scales subtended by the disk, but all momenta I and I of the 
angular-planar spectra are affected, with the net effect that they gain a positive/negative bias. 

This type of coherent behavior (an apparent bias) is in fact exactly what we observe in the following cases: (x 2 )2, 
(% 2 )g and, to a lesser extent, (x 2 )j ~ see Table n] The angular quadrupole i — 2, as well as the angular octopole 
I = 3, have all positive planar spectra (for all values of I which we were able to compute), indicated by probabilities 
larger than 50%. The planar hexadecupole I = 4 also has 8 out of 9 angular spectra assuming positive values (only 
I = 5 is negative). Again, we stress that the data analyzed in this Section relates to the full-sky maps, which are 
certainly still affected by residual galactic foregrounds, so it is important that we only draw conclusions after looking 
into the same results for the cut-sky (masked) maps - the subject of our next Section. 



B. Masked maps 



Our objective here, for completeness, is to perform the search for planar signatures in foreground-cleaned masked 
CMB maps. For this we will compare the masked (observational) variables y with the masked (theoretical) variables 
x. Following the WMAP science team recommendation to reject residual foregrounds contamination in temperature 
analyses of the CMB data we shall use the KQ85 mask, which cuts 18.3% of the sky data. Additionally, we also test 
the sensitivity of our results under the effect of the more severe KQ75 mask, which represents a cut-sky of 28.4%. 
Our final results for such analyses are shown in Tables (III I and (IV) and Figs. Q and 
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Regarding the results of our analysis, we start by noting that the sectors (l,£) — (2,5), (4,7) and (6,8), which 
have relative probabilities of 6.1%, 6.8% and 5% in the full-sky maps case, have now relative lower bounds given by 
23.1%, 20.9% and 11.1% for the mask KQ85 and 43.7%, 21.5%, and 14.6% for the mask KQ75, respectively. This 
reinforces an earlier claim that the planar anomalies detected in these scales are a consequence of residual foreground 
contaminations in the full-sky maps [56;. 

On the other hand, the angular scale I — 5 at I — 10 is still below the 10% threshold, even when we apply the 
severe KQ75 mask. Curiously, this scale was also reported by other groups |6j |66] as being anomalously spherical. 
The 10% level we found here corroborates this outcome, since a very negative value of (x 2 )g indicates a lack of planar 
modulations. 

The reader may have noticed that the sectors (10,8) and (2,9) in the KQ75 case are also below the 10% threshold. 
We argue that these results could be showing that the foreground cleaning procedures (used to generate the full-sky 
CMB map) may not only create false anomalies, but also hide others, in this case the sticky points at the angular 
scales I = 8, 9, which became clear only after the application of the recommended cut-sky mask. Nevertheless, whether 
these anomalies are real or caused by some systematic error in the map-making process, they suggest that further 
analyses are welcome. 
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76.1% 


61.9% 
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75.2% 
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46.2% 
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56.3% 
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16.9% 


88.4% 
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82.2% 


77.6% 


20.3% 


43.1% 
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47.0% 


37.1% 


29.7% 82.0% 


10 
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38.1% 


52.8%, 


3.2% 


60.8% 


60.4% 


36.7% 


32.5% 


68.9% 



Table III: Final probabilities of detecting, in a random ACDM universe with a galactic cut of 18.3% (KQ85 mask), a chi-square 
value smaller or equal to (x 2 )e(obs)- 
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Figure 4: Pdf's for (x 2 )e(th) (blue histograms), (x 2 )^(obs) (P ur P' e histograms) and for the difference (X 2 )f(th) — (x 2 )^(obs) ( s °lid 
red line) for maps with mask KQ85. Again, we only show some representative figures for simplicity. The final probabilities are 
shown in Table (III I. 



Finally, we turn to the coherence of the angular-planar spectra (x 2 )^ over a range of I's and ^'s. First, notice 
that the angular quadrupole £ = 2 shows again excessively positive modulations both for the KQ85 and KQ75 masks 
-although for the KQ75 case the planar / = 2 spectrum is marginally negative. The angular octopole £ = 3 is now 
only slightly anomalous: four out of the five spectra we computed turn out to be positive. The other case we discussed 
in the full-sky maps, I = 4, is still marginal: eight out of the nine spectra assume positive values. 

However, the most prominent feature that appears more clearly in the cut-sky maps is that £ — 5 has an anomalous 
coherence, with mostly rather negative values. In fact, £ — 7 and £ — 10 also show a high degree of uniformity with 
larger-than-expected modulations. As we showed in a previous Section, this type of coherence is the hallmark of a 
preferred plane (actually, of a disk-like region around a preferred plane.) Although we have not produced a test for 
this coherent signal, we believe our findings show that these angular scales deserve further investigation, and, given 
that they appear at approximately the angular size of the galactic cut, the coherence in (x 2 )^ could be pointing 
towards residual contaminations in the cut-sky maps. 
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38.3% 


52.2% 



Table IV: Final probabilities of detecting, in a random ACDM universe with a galactic cut of 28.4% (KQ75 mask), a chi-square 
value smaller or equal to (x 2 )e(obs)- 
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V. CONCLUSIONS 



We have used a top-down approach to search for deviations of statistical isotropy in the WMAP data by means of 
the recently introduced \ 2 test for angular-planar modulations [56 . 

These tools, designed to account for the presence of physical/astrophysical planes in temperature maps, can unveil 
any kind of planar modulations in these maps in a completely model-independent way. Furthermore, since improper 
subtraction of foreground contamination around the galactic plane in CMB data may lead to planar signatures, the 
angular-planar statistic can be a valuable tool to test the robustness of map cleaning procedures. Due to its generality 
and model independence, the angular-planar power spectrum C| m can be applied to any given temperature map -in 
fact, to any map whatsoever on the sphere S 2 . 

In this work we have applied this estimator to three classes of WMAP maps, namely, full-sky cleaned maps and 
masked cleaned maps due to the application of the WMAP 5-year KQ85 and KQ75 masks. We have included in our 
analyses an estimation of residual foreground contamination as well as an assessment of these types of errors in the 
final probabilities we calculated. 

The analysis with full-sky maps have shown that the angular quadrupole £ = 2 has a consistently positive planar 
modulation throughout the range I S [2, 10]. This systematic positivity could be an indication of a planar signature 
around £ = 2, although we would need to consider a wider planar range to confirm this suspicion. On the other 
hand, we have found slightly anomalous values in the angular-planar (l,£) sectors (2,5), (10,5), (4,7) and (6,8), with 
relative probabilities of 6.1%, 8.7%, 6.8% and 5%, respectively. 

As argued in [56], these anomalies are not drastically low, and can be attributed to residual foreground contami- 
nation. Our analyses with masked maps confirms that this is indeed the case for the sectors (l,£)=(2,5), (4,7) and 
(6,8). While in the full-sky map analysis the lower bound for the above probabilities was of ~5%, in the case with 
the mask KQ75 this bound was raised to 14.6%. This confirms earlier claims that foreground contamination may be 
an important issue when analyzing CMB anomalies, and that more reliable results should always be based in masked 
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maps analyses. 

On the other hand, we have found a significant trace of anomaly in the sector (I, £) = (10,5). This value seems to 
be robust under the effect of the two masks we used - lower values being allowed with only 4.6% of chance in the case 
of the mask KQ75. Curiously, the angular scale I = 5 was reported to be highly spherical by other groups [6j [66] . 
while we have found that it has a very low planar modulation. However, the multipole vectors used in that analysis 
do not represent physical directions in a straightforward manner, so the term "spherical" in this context should be 
used carefully. Nonetheless, it is interesting that this scale has conspicuously appeared as anomalous in two rather 
different tests of planarity. 

Perhaps the strongest indication of a possible anisotropy in the CMB maps was provided by the coherence of our 
X 2 test on the angular scales corresponding to I — 5, 7 and 10. As we showed, this uniformity of mostly negative 
or positive values of x 2 is equivalent to the signal left by a preferred plane (actually, a preferred disk with some 
thickness related to I), This may point towards contaminations in the masked maps which were left over from the 
foreground-cleaning process, although further work will be necessary to confirm this suspicion. 

As a final comment, we would like to mention that the approach we used to calculate our final probabilities, 
namely, the convolution of the observational and theoretical probability densities, is more rigorous than the widely 
used method of just computing the area under the theoretical probability density function. This approach generalizes 
the former, and implements in a natural way the uncertainty inherent in CMB observations, which in turn have a 
non-trivial impact in the final probabilities predicted by a specific cosmological model. 

In what regards prospects for further developments, we would like to point out that the tests used here are completely 
general and have a wide range of applications, among which we can mention polarization maps, tests of non-Gaussianity 
and even stacked maps of cosmic structure, like the galaxy cluster catalog 2Mass [67 . 
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